Entanglement of a qubit with a single oscillator mode 
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We solve a model of a qubit strongly coupled to a massive environmental oscillator mode where 
the qubit backaction is treated exactly. Using a Ginzburg-Landau formalism, we derive an effective 
action for this well known localization transition. An entangled state emerges as an instanton in 
the collective qubit-environment degree of freedom and the resulting model is shown to be formally 
equivalent to a Fluctuating Gap Model (FGM) of a disordered Peierls chain. Below the transition, 
spectral weight is transferred to an exponentially small energy scale leaving the qubit coherent but 
damped. Unlike the spin-boson model, coherent and effectively localized behaviors may coexist. 

PACS numbers: PACS numbers: 03.65.Yz, 71.27.+a 



The subject of quantum computing has led to renewed 
interest in the theory of quantum decoherence and also 
exposed a practical side to the theory. Scalable, persis- 
tent current designs for qubits involve measurement de- 
vices that are permanently coupled, leading to a contin- 
uous dephasing of the qubilii^. Before a quantum mea- 
surement can be made, the qubit must become entan- 
gled with the measurement device. The conflicting de- 
mands of quantum measurement and minimal dephasing 
is an active topic of study within the theory of mesoscopic 
quantum detector o 3 : 4 : 5 : 6 . 

In this letter, we study the problem of how a qubit 
becomes entangled with a single environmental mode 
when the qubit backaction on the environment is not 
neglected. This is perhaps the simplest quantum envi- 
ronment, but already displays great complexity. Models 
of dissipative quantum environments such as the spin- 
boson model (SBM) are intentionally weakly coupled and 
therefore neglect backaction. However, environments im- 
posed by quantum detectors may be strongly coupled 
and exhibit sharp spectral features. The entanglement 
described here might be realized in any environment (e.g 
detector or qubit-qubit coupler) that exhibits a low en- 
ergy monochromatic spectrum; as an example, we dis- 
cuss the environment imposed by an underdamped DC- 
SQUID detector. Lastly, we point out the equivalence 
between the present model and a particular model of dis- 
ordered fermions. 

Consider a qubit coupled to a single harmonic oscilla- 
tor described by the Hamiltonian 



H = Acr x + A 
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If the oscillator is replaced by a classical, adiabatic os- 
cillator (the limit ui — > 0, mu 2 — ► const), it is well 
knowni that JU exhibits a bifurcation at A 2 = Amw 2 
in that the energy of the ground state is now minimized 
by assuming a nonzero value of the oscillator displace- 
ment, x = ±xq- In the doubly degenerate ground state, 
the qubit is localized with a nonzero expectation value, 



(ct z ) 7^ 0. The classical calculation begs two questions: 
1) If the oscillator localizes the qubit, is all coherent be- 
havior destroyed? 2) If the qubit and oscillator weakly 
interact (A 2 << Amw 2 ) aren't they always entangled to 
some degree? 

We show that the onset of entanglement with the oscil- 
lator displacement becomes discontinuous in the massive 
limit, with a component of the ground state wavefunction 
playing the role of an order parameter in a second order 
phase transition. Below the transition, damped coherent 
behavior and entangled, effectively localized behavior co- 
exist. 

Entanglement of the qubit in our case is distinct from 
that of weak-coupling theories (such as the SBM) in that 
the ergodicity of the environment is effectively broken 
and fluctuations become non-Gaussian. This latter ef- 
fect, of course, cannot be seen by integrating out the en- 
vironmental degrees of freedom to find the effective action 
of the qubit. In this respect, our calculation bears some 
relation to the description of quantum measurement by 
Single Electron Transistor (SET) and overdamped DC- 
SQUID detectors 3 ^. In this scheme, the qubit degrees of 
freedom are integrated out to yield the response function 
of the detector. 

To represent the qubit, we choose the finite tempera- 
ture formalism of Popov and Fedotov 8 in which the action 
is written in terms of spinors satisfying modified bound- 
ary conditions. The imaginary time action corresponding 
to JTJ is 



S 



dTiP{T){d T + Aa x + \x(t)<t z )4>(t) + So (2) 



where S = J dr {^mx 1 if) + \kx 2 (r)) , k 
i/>(t) is a spinor 



U|(t) 



(3) 



The spin-1/2 is formally represented as a fermion 
spinor with an imaginary chemical potential. The latter 
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eliminates unphysical fermionic states from the Hilbert 
space. A gauge transformation removes the chemical 
potential and shifts the Matsubara frequencies to uj n = 
2n(n+ 1/4)//3. Conventional finite temperature field the- 
ory techniques may now be applied to the action (j2j; 8 - 
The fermions are now integrated out to yield an effective 
boson action: 



Z = I DipDipDxe~ s oc / Dxe 



X 



(4) 



S c fi = -trlog(-G + -^=x(w n - w m )cr z ) + S (5) 

where Go = (iui n — Aer^) -1 is the free fermion Green's 
function. The coordinate x(r) now represents a collective 
coordinate for a z and the oscillator displacement. 

Expanding the action to quartic order in A and re- 
grouping into dynamic and static terms, respectively, we 
obtain: 



+ ^E( fc +^°))i^«)i 2 ( 6 ) 
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where the coefficients are given by 
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Owing to the modified Matsubara frequencies, these 
sums are computed with a modified Fermi distribution, 
f(u>) = (ie^ + Defining the dimensionless control 

,2 

parameter a = the effective action is transformed 
back to imaginary time at zero temperature. 



Now we must find the corresponding dynamics of the 
qubit, in particular, the spin-spin correlation function. It 
may be shown that 

(T t (t 1 (t)ct j (0)) = i J DxiTr'ipai^ajif}} S{[x] e~ Sa{{ 

( 10 ) 

where S{[x] denotes averaging over the fermionic action 
for a particular realization of the boson field x(r). By 
Wick contractions, this average is related to the single 
fermion Green's function, in the presence of an inhomo- 
geneous field x(r), which solves the equation of motion: 
(d T + Aa x + Xx{t)<j z )G{tt';x(t)) = S(tt'). This ap- 
proach has been used to implement bosonization in 1-d 
fermion systems^. In the present case, the Green's func- 
tions must satisfy the boundary conditions implied by 
the shifted Matsubara frequencies: G(r + f3) — —iG(r). 
Furthermore, the boson field x must be written as the 
sum of an instanton trajectory and a small oscillation: 
x(t) = x\(t) + r(r). We now make the non-Abelian 
Schwinger ansata^, 



G(tt';x(t)) = U(r)g(TT';x i (T))U- 1 (T'), 



(11) 



where g satisfies the equation of motion but with the field 
x(t) restricted to the instanton trajectory: (d T + Aa x + 
Xxi(T)a z )g(TT' ; Xi(r)) — 5(tt'), and X[(t) — ±xq. Now 
U(t) must satisfy the auxiliary condition 



d T U+ [H,U] = ~Xr{r)a z U 



(12) 



where H = Aa x + Xxi(r)a z . 

A solution of the equation of motion for g which sat- 
isfies the boundary conditions is obtained for all r away 
from the instanton kinks (consistent with the dilute gas 
approximation); for r < t' 
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The mass is enhanced by a dynamical quantity and, most 
significantly, the action exhibits an instability at a — 
1. When a > 1 and u = 0, the action is minimized 
for x = ±A(2s^2)i/2^ For w 0; there is no true 
broken symmetry, however an approximate calculation 
of the instanton action 



Si = V2 
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shows that, although the symmetry changes at a > 1, the 
instanton is only stabilized for A/lo >> (a— 1)~ 3 / 2 . This 
defines the entangled regime of the qubit and oscillator 
in that the ground state is a coherent superposition of 
two states of the collective coordinate, x(t) = ±xq. 



where A is the Rabi energy A = \/A 2 + A 2 Xq. The 
Green's function g(r > t') is obtained by the replace- 
ment f(z) — > —f(z) where f(z) = /*(— z). The aux- 
illiary condition l|12|) leads to a set of Riccati equa- 
tions for an appropriate parameterization of U which 
are difficult to solve generally^. In the strong cou- 
pling limit (i.e. to lowest order in A/ A) the inter- 
acting Green's function is found to be G(tt' ; x(t)) = 

.g(rr';a; i (T))exp(-A(T z f T r(r)dr). 

Within this scheme we need to evaluate c zz (t — r') = 
(T t <j z {t)<j z (t')) s { [x]j other correlation functions are ob- 
tained similarly. c zz (t — r') involves two contractions 
yielding a DC part of c zz and a part at energy 2A, re- 
spectively. These correlation functions must in turn be 
averaged over the effective boson action, S e ff. 

(T T a z (T)a z (T')) = ( lim tr^G^r; x(ri)) 
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x lim trcr z G(T 2 T';x(T 2 )))s cff (13) 

+ (tra z G(rr'; x(t))<t z G(t't; x(t'))) So « 

The first contraction in l|13[) will involve products of 
Xi(r) at different times, giving the proper long time be- 
havior; the second contraction will involve the average 
over boson fluctuations r(r) that dress the qubit oscilla- 
tions at energy A. To evaluate the average over S e f{ we 
need both boson correlation function. Within the instan- 
ton approximation, the first contraction in (|13[l involves 



{TrX^Xiir')), 



-T\t-t'\ 



(14) 



where T — lo^J S , ;/27rexp —S\. The second contrac- 
tion in 1|13|) involves the average over fluctuations of 
r(r), (cosh J dTr(T))s C f f - Averages of this form and 
their fourier transforms have been discussed in several 
placesii. The relationship to the FGMi£ may now be 
seen: Replacing (in refiiS) the 1-d spatial degree of free- 
dom in the FGM by imaginary time, eqn. (J2J is the 
action for a fermion of frequency iA propagating in a 
Peierls chain with a gap function Xx(t). Averaging (in 
eqn. I|13|ll over the instanton fluctuations of x(t) with 
correlator Ijl4(l is equivalent to the disorder average in 
ref^fi with a correlation length given by the instanton 
time, r _1 . In the FGM, the density of states at the Fermi 
surface is accessed by the zero frequency limit uj — > i0 + 
which corresponds to A — * + , the strong coupling limit 
in the present work. The two models differ in that the 
correlation length of the gap disorder in the FGM is an in- 
dependently controlled parameter whereas the instanton 
fluctuations are generated spontaneously in the present 
model. In addition, the fermions in the present model 
satisfy twisted boundary conditions. 

Continuing to real frequencies and denoting the imag- 
inary part of the fourier transform of Ijl3(l by s Z z{v), we 
obtain for a > 1 a correlation function of the form 



+ ^-(x )2 + ( f )2) 



(15) 



n n+m — 2p 
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for v > where v = a—. This is our main result. 
Equation l|15[) shows that for the broken symmetry phase 
(a > 1) two distinct spectral features are formed. The 
low energy feature with weight u> cnt = ^(Xxq/ A) 2 ~ 
^— - (at energy T) corresponds to the qubit entangled 
with the oscillator and in a superposition of the form 
| + ; +xo) ± |— ; — xo), where | + / — ) refer to states of the 
qubit with nonvanishing z-polarization. The high en- 
ergy set of delta functions (at energy ~ 2 A) , corresponds 
to the decoupled (unentangled) but dephased qubit with 
weight iffroe — 1/a. Although the spectrum is discrete, 
dephasing may be estimated by considering the envelope 
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FIG. 1: Numerical diagonalization of Hamiltonian Q. Dy- 
namical spin-spin correlation function showing the emergence 
of the low energy instanton feature corresponding to entangle- 
ment. In this calculation, ui = 0.01, A = 0.5, k = 0.2. Peaks 
are artificially broadened and vertical scale on curves for the 
largest four values of a is magnified x2. 



of the sidebands associated with the primary resonance. 
For small u> (large p) the Poisson distributions become 
approximately Gaussian and the dephasing rate is es- 
timated to be up 1 ' 2 ; the fractional width of the reso- 
nance, 1/Q, is then (a-^) 1 / 2 . The qubit becomes criti- 
cally damped when this factor is unity. When m — > oo, 
keeping k and therefore a constant, the width of the res- 
onance goes to zero as 0(m -1 / 4 ) (although the number 
of bosons diverges) , in agreement with the classical cal- 
culation. 

The spectral weights appearing in equation l|15f) sat- 



isfy the sum rule: 
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1. At critical coupling, 



a redistribution of spectral weight is initiated and spec- 
tral weight flows from the coherent feature to the entan- 
gled one. Unlike the SBM where the system spin be- 
comes overdamped before it becomes localized, coherent 
and effectively localized behavior coexist in the present 
mode!—. An approximate numerical diagonalization of 
the Hamiltonian Q confirm these results. Figure 1 shows 
that a feature of energy T appears in the correlation 
function s zz (y) for a ~ 1, while the primary resonance 
remains distinct from the instanton feature. As expected 

1/2 

from the effective action S c g, w c ^ t is the order param- 
eter for a second order phase transition; in the strict 

1 /2 

adiabatic limit (A/w — > oo), the derivative dw c ^ t /da 
would diverge. To illustrate the sharp onset of entan- 
glement, we have computed the entanglement entropy 
S c = — tr / olog 2 p, where p is the ground state reduced 
density matrix of the qubit. The entropy, shown in fig. 
2, is seen to sharply increase from zero at the onset of 
entanglement (a = 1). Similarly to the order parameter, 
the slope of this graph appears to have a discontinuity 
at a = 1 as A/oj — > oo. This is not unexpected, not- 
ing that thermodynamic entropy is given in terms of the 
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FIG. 2: Entanglement entropy computed numerically from 
Hamiltonian with w = 0.002, A = 0.5, k = 0.2. Solid line 
is a plot of Xq scaled by a constant factor. 

Gibbs free energy by 5th = —dG/dT. Taking G to be 
the static limit of 5 e ff, —dG/da — \kx\, which is pro- 
portional to the square of the order parameter. Thus, 
entanglement entropy appears to behave analogously to 
thermodynamic entropy in this second order phase tran- 
sition. 

As a potential example of this transition, we now 
turn to experimental realizations of superconducting 
qubits. Ini it was demonstrated that the dominant 
source of decoherence in a persistent current (phase) 
qubit was the electromagnetic environment of the DC- 
SQUID, characterized^ by an RLjC impedance Z(oj) (Lj 
is the Josephson equivalent inductance.) Dephasing fol- 
lows from the spin-boson model where J(u>) oc — ReZ(w); 
in particular, J is ohmic at low frequencies. However, 
J (to) is highly structured, dominated by the plasma mode 
at the frequency u> = 1/s/LjC with a Q ~ 10 for typ- 
ical parameters^. The a z coupling to the qubit comes 
from the circulating component of current in the DC- 
SQUID. The effective coupling constant, A, reflecting 
the mutual inductance (M) between the persistent cur- 
rent in the qubit (I p o~ z ) and circulating current of the 
DC-SQUID is then A = 2MI I P sin ±^sin#, where I is 
the critical current, is the net Josephson phase and 
4> = 27r($/$o) the dimensionless flux at the operating 
point. Linearizing the effective potential of the SQUID 



about the operating point, the stiffness is found to be 
k = 2& J 4/q cos 2 \<i> — if, where I e is the external cur- 
rent bias. Thus the effective potential softens and the 

A 2 

dimensionless coupling constant a = 4^ increases as I e 
approaches the switching value. 

To illustrate our point, we use M = 20 pH, I p = 900 
nA, Jo = 150 nA, and fix I e at 0.95 of the critical switch- 
ing value, 2Iq cos ^</>, to obtain a dimensionless coupling 
of a ~ 1. Once entangled, two distinct coherent oscil- 
lations in a z would be observed: the primary oscillation 
at 2 A and the entangled one at V. Within this model, 
the oscillation at T is undamped; however, external cou- 
plings would be expected to cause dephasing. The single 
mode contribution to the dephasing of the 2A oscilla- 
tion becomes significant when p(= a^j) >> 1 in which 
case Q = (^j) 1 ^ 2 - Using a — 1 and the adiabatic ra- 
tio A/ui = 10, yields a Q for the qubit of Q ~ 3.3. As 
in fig. 1, however, the spectrum is discrete consisting 
of several sidebands. In contrast, the ohmic spectrum, 
J(lo) = rju>, for the same DC-SQUID parameters (and 
lead resistance R = 10017) has a dimensionless coupling 
of r] ~ 0.03. The ohmic contribution to the dephas- 
ing rate is 1/r ~ ^ 7T7 1 i": which vanishes for both 
zero bias and zero temperature (the present case). For 
e = A ~ 1 GHz, the ohmic Q becomes comparable for 
a temperature of T ~ 10 mK and diverges for T = 0. 
In contrast, the single mode contribution to dephasing 
remains constant at T = 0. 

In conclusion, we have considered a model in which a 
qubit is strongly coupled to a single environmental mode 
and described the resulting dephasing and entanglement. 
We have solved the resulting model by appealing to a for- 
mal analogy with the Fluctuating Gap model, a model for 
disordered fermions. We applied these results to an un- 
derdamped SQUID environment, however, this transition 
might be studied in a variety of existing qubit designs, 
the main attribute being a variable coupling between the 
qubit and an adiabatic (A/ui ~ 10) oscillator mode. 

We thank Jeffrey Clayhold, Boris Fine, Philip Stamp, 
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